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We present a preliminary set of updated NLO parton distributions. For the first time 
we have a quantitative extraction of the strange quark and antiquark distributions and 
their uncertainties determined from CCFR and NuTeV dimuon cross sections. Addi- 
tional jet data from HERA and the Tevatron improve our gluon extraction. Lepton 
asymmetry data and neutrino structure functions improve the flavour separation, par- 
ticularly constraining the down quark valence distribution. 

There are many reasons why an update p] to the MRST2004 parton distributions [2] 
is necessary. The MRST (now MSTW) group have used a general-mass variable flavour 
number scheme (VFNS) since 1998 [3], but a new scheme [4], although primarily devised for 
use at NNLO [5], also changes the details of the NLO implementation. We now make use 
of the f astNLO package [BJ, which provides an efficient combination of NLO partonic cross 
sections with parton distributions, allowing for a more rigorous inclusion of jet data into 
global parton fits. 

There are also many new data sets avail- 
able to include: CHORUS and NuTeV neu- 
trino structure functions Q 2 ) and 
F^{x,Q 2 ) [7], CCFR and NuTeV dimuon 
cross sections [8] providing a direct con- 
straint on s and s, CDF Run II lepton asym- 
metry data [9] in two different Et bins, 
HERA jet data [TU], and new CDF Run II 
high-Er jet data [TT]. We also include all 
recent charm and bottom HERA structure 
function data [12]. 

Let us consider the effects of some of 
these new data sets. It has already been 
noticed that the NuTeV neutrino struc- 
ture functions are not completely compat- 
ible with the older CCFR data, both by the 

experiment themselves and by groups performing fits [131 114] . Previous parton distributions 
were perfectly compatible with the CCFR data when using EMC-inspircd <5 2 -indepcndcnt 
nuclear correction factors, whereas the NuTeV structure function data are difficult to fit at 
high x. Moreover, recent CHORUS structure function data (which were taken with a lead 
rather than an iron target) turn out to be fairly consistent with the CCFR data. Thus, there 
is a problem at high x. However, the partons in this region are already very well determined 
from charged lepton structure functions. The important information from charged current 
data is in the region x < 0.3, probing the separation into valence quarks and sea quarks, and 
there is general consistency between data sets here. Hence, we choose to exclude from the 
fit all neutrino structure function data for x > 0.5. In our analysis we now also implement 




Figure 1: The distribution xd v (x, lOGeV 2 ). 
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Figure 2: The input distributions s + (x,Qq) (left) and s (x,Qq) (right). 



more sophisticated flavour-dependent nuclear correction factors [TS] extracted using NLO 
partons, allowing a ~ 3% uncertainty on these factors. 

The CDF lepton asymmetry data constrain similar partons as the neutrino structure 
functions. The W asymmetry at the Tevatron is defined by Aw{yw) oc da(W + )/dyw — 
da(W~)/dyw ~ u{x\)d{x2) — d{xi)u{x2) , where Xi t 2 = (Mw/Vs) exp(±yw)- In practice 
it is the final state leptons from the W decay that are detected, so it is really the lepton 
asymmetry A(yi) = [a(l + ) — a(l~)]/[a(l + ) + (r(l~)] which is measured. Defining the emission 
angle of the lepton relative to the proton beam in the W rest frame by cos 2 9* = 
leads to yi = yw ± l/21n[(l + cos 0*)/(l — cos 8*)]. To a good approximation 



a(l + ) -a(l ) oc u(xi)d(x 2 )(l - cos 9*) 2 + d(x 1 )u{x 2 ){l + 
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where the antiquark term is boosted by the large (1 + cos 9*) 2 and the asymmetry is fairly 
sensitive to antiquarks at lower Et- Hence, it probes the separation into valence and sea 
quarks, particularly for the less well-constrained down quark. We find that the CDF Run II 
data influence d(x,Q 2 ), and that there is some tension with the neutrino structure function 
data. Our fit gives x 2 = 41 for the 22 lepton asymmetry data points, although much of 
this x 2 comes from only 2 data points. Overall, dy(x, Q 2 ) has a slightly different shape to 
previously [16j . seen in Fig. [TJ The uncertainty is generally bigger than before, despite more 
constraints, due to a better parametcrisation when determining uncertainty eigenvectors. 

The CCFR and NuTeV dimuon cross sections for neutrino (antineutrino) scattering, 
v^{p^)N — > arc sensitive to the strange quark (antiquark) distribution through the 

LO partonic process v^s{v^s) — > c[i~ (c/j, + )X . In previous MRST fits, CCFR dimuon data 
have been indirectly used to justify the input parametcrisation assumption of 

s(x, Ql) = s(x, Ql) = ^[u(x, Ql) + d(x, Ql)\ (k r» 0.5), 



at the input scale of Qq = 1 GcV 2 . We now relax this assumption and fit the CCFR/NuTcV 



dimuon cross sections directly. Defining s + i 
taken to be s + (x, Ql) = A + (1 - x) v + S(x, * 



+ s and s = s ■ 
and s~(x, Qq) 



s, then the input forms are 



-l+<5- 



(1 



x/xq), where S(x,Ql) is the total sea quark distribution and xq is determined by zero 
strangeness of the proton, i.e. L dx s~(x, Qq) = 0. The preference for extra freedom in both 



DIS2007 



s + and s~ is confirmed by the fit, with the global \ 2 improving by 15 when s + is allowed 
to go free, and by a further 14 when s~ is allowed to go free. There is no improvement with 
the addition of further free parameters. Indeed, A_ and 5— are highly correlated so we fix 
5- at 0.2, i.e. a valence-like value. 



CDF Run II inclusive jet data, % 2 = 56/76 pts. 
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The input s + and s~ distribu- 
tions from the best fit, with approx- 
imately 90% confidence level uncer- 
tainty bands, are shown in Fig. [2] 
There is a reduced ratio of strange 
to non-strange sea distributions com- 
pared to our previous default k = 0.5, 
and a marked suppression at high x, 
i.e. low W 2 , probably due to the effect 
of the strange quark mass. The inte- 
grated momentum asymmetry is pos- 
itive, 0.0028t8;°°i| at Ql = lGeV 2 , 
decreasing to 0.0021±°;°°™ at Q 2 = 
10 GeV 2 . The results on both s+ and 
s~ are qualitatively consistent with 
those obtained by the CTEQ group 
[17] . Directly fitting the s and s dis- 
tributions affects the uncertainties on 
the light quarks. Until recently all par- 
ton fitting groups assumed s + to be a 
fixed proportion of the total sea in the 
global fit. An independent uncertainty 
on s and s feeds into that on the u 

and d quarks, since the neutral current DIS data on F2(x,Q 2 ) constrain the combination 
4/9(u + u) + l/9(d + d + s + s). Consequently, the size of the uncertainty on the sea quarks 
at x - 0.001 - 0.01 at hard scales Q 2 - M 2 ? roughly doubles from - 1.5% to - 3%. 

We now use the fastNLO [B] package 
during the fit to calculate Tevatron and 
HERA jet cross sections at NLO, improv- 
ing upon previous approximations for the 
former and allowing our first analysis of 
the latter. We find a slight change in the 
shape of the gluon distribution using only 
Tevatron Run I data (for which we now in- 
clude previously absent hadronisation cor- 
rections). The fit to the HERA jet data 
is excellent, although we do not find that 
they provide a strong constraint within the 
global fit. We also include CDF Run II in- 
clusive jet data jTTj obtained in different ra- 
pidity bins using the kr jet algorithm. We 



Figure 3: The fit to the CDF Run II jet data. 
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Figure 4: The distribution xg(x, 2 GeV 2 



obtain a very good fit, \ 2 = 56 for 76 data points, as seen in Fig. [3l at the expense of a 
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slight deterioration in the fit to the D0 Run I inclusive jet data [TS]. The CDF Run II data 
prefer a smaller gluon distribution at very high x compared to the Run I data included in 
the MRST2004 fit [2]. The gluon distribution is shown in Fig. g] compared to MRST2001 
|16j . The uncertainty is similar at high x, but is actually larger overall due to allowing an 
additional parameter to go free when determining the eigenvectors. The uncertainty on the 
gluon is extremely large at x = 10 -5 due to our lack of theoretical prejudice on the form 
taken in this region. 

Finally, we investigate the dependence of 



the fit on m c , the pole mass of the charm 
quark. We take the pole mass of the bottom 
quark to be m& = 4.75 GeV. We vary m c 
in steps of 0.1 GeV, finding a preference for 
m c = 1.4 GeV with an uncertainty ~ 0.2 GeV. 
The results are shown in Table [TJ There is 
a clear correlation between m c and a,g(M§). 
The value of the (NLO, MS) coupling for the 
best fit is as(Mz) = 0.1206, very similar to 
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that in MRST2004 g]. Ta bl c 1: The dependence of the fit on m c . 

Overall, the inclusion of new data and the 
changes in the analysis have a significant impact on our NLO parton distributions. Based 
on the developments in this paper we will soon have fully updated NLO and NNLO partons 
for the LHC complete with uncertainties. 
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